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Abstract 

We present results on measurements of jet shapes and jet azimuthal decorrelation 
from pp collisions at = 1.8 TeV using data collected during the 1992-1993 run 
of the Fermilab Tevatron. Jets are seen to narrow both with increasing Et and 
increasing rapidity. While HERWIG, a parton shower Monte Carlo, predicts slightly 
narrower jets, it describes the trend of the data well; NLO QCD describes qualitative 
features of the data but is sensitive to both renormalization scale and jet definitions. 
Jet azimuthal decorrelation has been measured out to five units of pseudorapidity. 
While next-to-leading order QCD and a leading-log approximation based on BFKL 
resummation fail to reproduce the effect, HERWIG describes the data well. 



^Presented at Les Rencontres de Physique de la Vallee d'Aoste: Results and Perspectives in Particle 
Physics, La Thuilc, Italy, 3-9 March 1996. 



1 Introduction 



As progress has extended perturbative quantum chromo dynamics (QCD) beyond lead- 
ing order, theoretical comparisons of a more detailed nature are being performed with 
experimental data. The advent of both next-to-leading order (NLO) and leading-log 
approximations (LLA) has enabled predictions for both event topologies involving more 
than two jets and the internal structure within a jet. This paper reports on two such stud- 
ies, one involving a detailed examination of the transverse energy flow within jets and the 
other involving the first study of azimuthal decorrelation of jets with large separations in 
pseudorapidity. 

In both studies, NLO QCD is the first order in which a theoretical prediction is 
possible. Although NLO predictions have been successful in describing inclusive single jet 
and dijet cross sections, deviations from the measured distributions or a large sensitivity 
to the choice of renormalization scale may indicate large higher order corrections. In this 
case, an all-orders approximation in the form of a LLA may produce better results. For 
the jet shape analysis, a model where only one additional radiated parton is allowed 
may prove inadequate to reproduce all aspects of the transverse energy flow. Additional 
fragmentation effects may indeed be important. In jet production with large rapidity 
separations, the presence of multiple scales in the calculation may require a resummation. 
Two alternative methods are explored, one resumming InQ^ via DGLAP splitting 
functions and the other resumming ln{s/Q'^) using the BFKL ||^ equation. A comparison 
of experimental results to theoretical predictions will determine the necessity of including 
such higher order corrections. 

2 Jet Shape 

The data for this study |^ were collected using the D0 detector which is described in 
detail elsewhere |^ . The critical components for this study are the uranium-liquid argon 
calorimeters. These calorimeters provide hermetic coverage over most of the solid angle 
with a transverse segmentation of 0.1 x 0.1 in t] x (j) [rj = — ln[tan(6'/2)], where 6 is the 
polar angle with respect to the proton beam and is the azimuthal angle). The measured 
resolutions are 15%/\/E and 50%/ (E in GeV) for electromagnetic and single hadronic 
showers respectively. 

The data were collected using four separate triggers, each requiring a minimum Ej' in 
a specified number of trigger towers (0.2 x 0.2 in 77 x 0) at the hardware level. Events 
satisfying these requirements were sent to a processor farm where a fast version of the 
offline jet algorithm searched for jet candidates above some Erp threshold. These events 
were used to populate four non-overlapping jet Ft ranges: 45-70, 70-105, 105-140, and 
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> 140 GeV. Each trigger was used only where it was fully efficient. 

Jets were reconstructed offline using an iterative fixed cone algorithm with a ra- 
dius R = 1.0 in 7] — space. The final jet angles were defined as follows: 6'jet = 
tan^^ ^ {T^iE^iY + {T^iEyiY) /lliE^i and 0jct = iim~^{T.iEyi/'LiE^i) with z in the beam 
direction and the sums extending over calorimeter towers whose centers were contained 
within the jet cone. These definitions differ from the Snowmass algorithm ?7jet = 
T,iETir]i/T,iETi and 0jet = ^iETi4>i/'^iETi- In both cases, the Et of the jet was defined 
as T^iETi- After a prehminary set of jets was found, overlapping jets were redefined. Two 
jets were merged into one jet if more than 50% of the Et of the jet with the smaller 
Et was contained in the overlap region. The direction of the new jet was defined as the 
vector sum of the two original jet momenta, and the energy was recalculated. If less than 
50% of the Et was contained in the overlap region, the jets were split into two distinct 
jets. In this case, the energy of each calorimeter cell in the overlap region was assigned to 
the nearest jet and the jet directions were recalculated. Two different rj regions {\ri\ < 0.2 
and 2.5 < \r]\ < 3.0) were studied. 

For this study we define the quantity p(r) = {1/ Njf,ts)^}etsET{r) / Et^t = 1) (A^jets is 
the number of jets in the sample) which is the average fraction of the jet transverse energy 
contained in a subcone of radius r. For a given value of r, a larger value of p represents 
a narrower jet. Individual subcones were corrected for underlying event and noise effects 
only. 

In Fig. I, the average Et fraction (p) versus the subcone radius is shown for the two 
different t] regions for jets between 45 and 70 GeV Et- Also shown in this figure are the 
predictions from HERWIG , a Monte Carlo prediction which incorporates higher order 
effects through a parton shower model. Detector effects were modeled using a detailed 
detector simulation based on GEANT p. Although properly describing the qualitative 
behavior of the data, in both the central and forward regions the theoretical predictions 
are narrower than the measured jet shapes. In order to remove detector effects, correction 
factors were determined using three different event generators, comparing the jet shape 
before and after detector simulation. In the subsequent distributions, the detector effects 
have been removed using these correction factors. 

In Fig. H, the jet shape for the four different Et regions are shown for central jets 
{\ri\ < 0.2). We see that as the Et of the jet increases, a larger fraction of its total 
Et is contained in the inner subcones, i.e. the jets are narrower. This result is in good 
agreement with CDF ^ measurements (Fig. ^ where jet shapes were determined using 
charged particle tracks only. In Fig. ^, the jet shape for a single Et range in both the 
central and the forward r] region is shown. We see that jets of the same Et become 
narrower as the jets become more forward. 

Also shown in Figs. |^ and |^ is a comparison to NLO QCD predictions made with 
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Figure 1: The average integrated Et fraction versus the subcone radius is plotted for 
the data and HERWIG Monte Carlo, before calorimetric effects are removed, for the Et 
range 45-70 GeV for (a) I77I < 0.2 and (b) 2.5 < I77I < 3.0. 
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Figure 2: The measured jet shapes, with calorimetric effects removed, compared to NLO 
predictions (JETRAD) with three rcnormalization scales for \ri\ < 0.2 for the jet Et range 
(a) 45-70 GeV, (b) 70-105 GeV, (c) 105-140 GeV, (d) greater than 140 GeV. 
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Figure 3: A comparison of the D0 and CDF jet shapes. The CDF measurement was 
made using charged tracks only for 0.1 < \r]\ < 0.7 and average Et = 100 GeV. The two 
D0 measurements (|?7| < 0.2), with detector effects removed, were made at E't's above 
and below that of CDF. Good agreement is seen. 




Figure 4: Comparisons of jet shapes from data and NLO predictions using different jet 
direction definitions for jets with 70 < Et < 105 and for (a) \ri\ < 0.2 and (b) 2.5 < |?7| < 
3.0. 
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JETRAD |jTo| using CTEQ2M [|Tl|] parton distribution functions (pdf ) and three different 
renormalization scales(/i): 2Et, Et, and Et/2. Studies sliowed the theoretical predictions 
were insensitive to choices in pdf. Since this is a NLO QCD prediction containing at most 
three partons in the final state, jet shape arises from combining two final state partons 
into one jet. For this study, two partons were considered to constitute one jet if their 
distance in rj - (p space from their vector sum was less than one. Using the D0 definitions 
of jet angles, central jets are found to be generally wider than predicted by JETRAD. 
Forward jets are narrower than predicted and the observed narrowing of forward jets with 
increasing Et is not reproduced by JETRAD. 

As shown in Fig. ^ using the Snowmass definition of jet angles generally improves the 
agreement between data and predictions. While the experimental jet shape distribution 
is only moderately affected by the choice of jet angle definition, the theoretical prediction 
in the forward region is very sensitive to this choice. Qualitatively, the NLO calculations 
using the Snowmass algorithm describe the measured jet shapes reasonably well. However, 
these predictions, being the first order in which a nontrivial jet shape emerges, show large 
sensitivity to the choice of renormalization scale with no single choice able to describe the 
data at all E^s. 



3 Jet Azimuthal Decorrelation 

In this study ||12| we are interested in jets widely separated in pseudorapidity and at 
relatively low Et- By measuring the correlation of the dijet system we gain information 
about the radiation pattern in the event. The data used for this study were collected in a 
similar fashion to the jet shape study. In this case, a single trigger was used. It required 
Et > 7 GeV at the hardware level and searched for jet candidates with an Et threshold 
of 30 GeV in the processor farm. Offline, jets were reconstructed using an iterative fixed 
cone algorithm with a radius of R = 0.7 in r/ — space. Both the final jet angles and the 
merging and splitting criteria were defined as described above. 

From the sample of jets collected with Et > 20 GeV and \ri\ < 3.0, the two jets at 
the extremes of pseudorapidity were selected. To remove any trigger inefficiencies, one of 
these two jets was required to have Et > 50 GeV. The relevant quantities used in this 
analysis are Arj = \rii — 772! and = 0i — 02 (where the subscripts refer to the two jets 
described above). 

Since the trigger was implemented out to \ri\ < 3.0, the maximum value of At] is 
~ 6.0. This is illustrated in Fig. ^ where we plot the pseudorapidity interval Arj for 
events passing our cuts. In Fig. ^, the azimuthal angular separation, 1 — Acp/ir is plotted 
for three selected unit bins of Arj centered at At] = 1,3, and 5. The correlation of the two 
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Figure 5: The pseudorapidity inter- 
val, At] = |?7i — 772!, of the two jets at 
the extremes of pseudorapidity. The 
coverage extends to Ar] ~ 6. The 
errors are statistical only. 



Figure 6: The azimuthal angle differ- 
ence, A0 = 01 — 02, distribution of 
the two jets at the extremes of pseu- 
dorapidity plotted as 1 — Acp/n for 
At] = 1, 3, and 5. The errors are 
statistical only. 



jets at the extremes in pseudorapidity is evident in the shape of this distribution which 
becomes wider (less correlated) at Arj increases. A correlation variable, (cos(7r — A0)), 
has been defined to study the effect quantitatively. This quantity is plotted in Fig. ^ as 
a function of Ar]. For the data, the error bars represent the statistical and uncorrelated 
systematic errors added in quadrature. An error band at the bottom of the plot represents 
the correlated systematic error due to the energy scale and detector effects. Also shown 
in Fig. 1^ are the theoretical predictions from HERWIG, NLO QCD as implemented in 
JETRAD, and BFKL resummation The errors for the predictions are statistical 



only. 

For the data, many sources of systematic error were investigated. For all values of 
At] the energy scale error dominates, except for Arj = 5, where the measurement is 
statistically limited. Other systematic errors due to the out-of-cone showering correction, 
angular biases present in the jet reconstruction, and our selection cuts are included. 
Although no attempt was made to correct back for detector effects, the size of any such 
effects was estimated using HERWIG events that were subjected to a detailed simulation 
of our detector based on GEANT. The size of the effect was negligible for Arj < 3 and 
only ~ 0.03 at Arj = 5 and are included in the systematic correlated error band. 

Systematic effects were also studied in the theory. The renormalization and factor- 
ization scale /X was varied in the NLO prediction from n = pip"-^ /2 to 2p^°'^. The largest 
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Figure 7: The correlation variable used in this analysis, (cos(7r — A0)) vs. Arj, for the 
data, JETRAD, HERWIG, and the BFKL calculations of Del Duca and Schmidt. 

variation seen in (cos(7r — A0)) was less than 0.026. Different parton distribution func- 



tions were used (CTEQ2M, MRSD- [y], and GRV [T|]) resuhing in variations of less 
than 0.003. In addition, the Monte Carlo study was repeated using the Snowmass jet 
angle definitions. The difference in (cos(7r — A0)) was less than 0.013. As a cross check, 
a subset of the data was analyzed using the Snowmass jet angle definitions as well. The 
resulting variations were less than 0.002. 

The data in Fig. ^ show a nearly linear decorrelation effect with pseudorapidity. For 
small At] both HERWIG and JETRAD describe the data well. However, at large Arj 
JETRAD, which is leading order in describing any decorrelation effect, shows deviations 
from the measurement predicting too little decorrelation. The BFKL predictions, valid 
for large agArj, are shown for At] > 2. As the pseudorapidity interval increases, this 
leading-log approximation predicts too much decorrelation. The HERWIG prediction, 
which includes higher order processes in the form of a parton shower model including 
angular ordering, describes the data well over the entire Arj region. 
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4 Conclusion 



We have presented results on both jet shapes and azimuthal decorrelation in dijet sys- 
tems. These resuhs have been compared to various theoretical predictions. We see that 
jets become narrower as their Et increases for a fixed r] and that they become narrower 
as their rj increases for fixed Et- The measured jet shapes are insensitive to jet angle 
definitions while the NLO QCD predictions are sensitive to jet angle definition. Reason- 
able behavior is seen when Snowmass jet angle definitions are used. We see that the dijet 
systems decorrelate in a linear fashion out to A77 = 5. Comparisons to various theoretical 
predictions show that NLO QCD predicts too httle while BFKL resummation predicts 
too much decorrelation. HERWIG reproduces the data over the entire A77 range explored. 
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